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Abstract. We consider gradient descent equations for energy functionals of 
the type S(u) = ^(u(x), A(x)u(x)) L 2 + Jq V(x, u) dx, where A is a uniformly 
elliptic operator of order 2, with smooth coefficients. The gradient descent 
equation for such a functional depends on the metric under consideration. 

We consider the steepest descent equation for S where the gradient is an 
element of the Sobolev space Hi 3 , (3 £ (0, 1), with a metric that depends on A 
and a positive number 7 > sup | V22 1 ■ We prove a weak comparison principle 
for such a gradient flow. 

We extend our methods to the case where A is a fractional power of an 
elliptic operator. 

We provide an application to the Aubry-Mather theory for partial differen- 
tial equations and pseudo-differential equations by finding plane-like minimiz- 
ers of the energy functional. 



1. Introduction 

In this paper wc prove a comparison principle for steepest descent equations, 
in the Sobolev gradient direction (see ©)■ When one is interested in minimizing 
functionals of the type 

S(u) = -{u(x),A(x)u(x)) L 2 + [ V(x,u)dx, (1) 
z Jn 

where A is an elliptic operator, it is natural to consider the gradient descent equation 
dtu = — \7S(u). The gradient of S depends on the metric under consideration. The 
metrics and Sobolev gradients we consider are explained in Sections 11.21 and 11.41 
See also }Neu97j . 

We are interested in comparison principles for gradient descent equations. One 
possible application is the use of comparison in establishing a priori bounds on the 
critical points of S in the context of the Aubry-Mather theory for partial differen- 
tial equations and pseudo-differential equations (see |dlLV09j ). but we develop the 
general methods in detail apart from applications. This will be discussed further 
in Section [5] 

More concretely, we consider a self-adjoint, uniformly elliptic operator A given 
by the formula 

d 

Au = - d *i (a i3 (x)d Xi u) = -div(o(s)Vu), (2) 
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where the coefficient functions, a y <G C°°(R d ) are symmetric in i,j and we have 
positive constants Ai, A2 such that for everyx £ K d 

d 

Ai|£| 2 < J2 aij ( x Mi < A ^l 2 > e M d . 

Then, for a suitably large constant 7, and for [3 £ (0, 1), we will show a comparison 
principle for the flow defined by the evolution equation 

d t u= -{ 1 + A) 1 -Pu+{ 1 + A)-e{ 1 u-V 2 {x 1 u)) 1 (3) 

where V(x, y) £ C r (R d x R), r > 2, and V2 denotes the derivative of V with 
respect to its last argument. The fractional powers of 7 + A that appear in Q 
will be defined in Section 13.21 As we shall show over the next few sections of this 
introduction, equation (J3J is the steepest descent equation for S in the Sobolev 
space H@ with inner product (u,v) H p = ((7 + A) l3 u,v) ^2 , as explained in Section 
If .21 The domain and boundary considerations for equation ([3]) are discussed in 
the following section. A sufficient lower bound for the constant 7 will be given in 
Section [4] and will depend on the nonlinear term V. 

1.1. Boundary Conditions. Our main results apply to Dirichlet boundary con- 
ditions for domains C R d that are compact with smooth boundary, as well as for 
periodic boundary conditions (i.e. £1 = T d = R d /Z d ). 

In the application in Section [5] wc will also consider Ct = NT d = R d /NZ d 1 for 
which the reasoning regarding N = 1 applies. More succinctly, we pose ([3]) as 
an initial-boundary value problem, with u(0,x) = uq(x) £ L°°(0) and one of the 
following two cases: 

u = on dfl (4) 

u(x + e,t)=u(x,t) Me elf 1 (5) 

In the periodic setting we will require the functions a lJ (x) and V(x,y) to have 
period one in all variables. 

Since our results rely on arguments that apply to both the Dirichlet and periodic 
settings, we will not distinguish between the different cases when stating the results. 

1.2. Sobolev Spaces. There are several equivalent definitions of the Sobolev spaces 
H^n), s e R. For the integer case, we take H m (tt) = {u € L 2 (tt) : D a u g 
L 2 (f2), V|a| < to}, then the intermediate spaces may be defined by interpolation 
methods |Lun95j , |Tay96| . Alternatively, one can use the Fourier transform to de- 
fine H s (R d ) = {u G L 2 (R d ) : (1 + \t\ 2 ) s/2 u e L 2 (IR d )} and then H s (n) = {u\ n : 
u E H s (M. d )}. The case where = T d is handled simply by replacing the Fourier 
transform with the Fourier series. The factor (1 + |£| 2 ) s / 2 makes it clear that the 
operator I — A is the foundation of these spaces. Indeed, u £ H s (M. d ) if and only 
if (/ — A) s / 2 u £ L 2 (M. d ), where (/ — A) -s / 2 is a particular case of the general defi- 
nition of the power of an elliptic operator given in equation |6]), and (/ — A) s / 2 is 
the inverse of (I — A) -5 / 2 . 

In fact, because 7 + A is an (order 2) elliptic, self-adjoint operator, we can define 
the Sobolev space H* A in the same manner as above, but replacing I— A with 7+^4, 
for 7 > 0. The inner product on A is given by (u, i>) s , 7 .a = ((7 + A) s u, v)o = 
((-f+A) s u 7 v) L 2. In [ShuOlj . page 57 it is shown that the topology on A generated 
by the norm obtained from the above inner product is identical to the standard 
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topology on H s (i.e. the topology generated by (u,v)h* = ((I — A) s u,v) L 2). Thus, 
we henceforth omit the subscripts 7 and A when referring to iff A and (•, -) s .-y,A- 
Will will write (•, -) s for the inner product on H s and || • || s for the norm on H s . 
Though the topologies are equivalent, the gradient of S depends on the chosen inner 
product. Thus, the gradient flow and therefore the comparison principle depend on 
the choice of inner product. 

1.3. Fractional Powers of Elliptic Operators. For s > 0, the operator (7 + 
A)~ s is self-adjoint, bounded, linear, invertible from H r to H r+2s , and is defined 

by 

(7 + A)- s = z- s ( 7 + A - z)- 1 dz, (6) 

2m J r 

where T is a rectifiable curve in the resolvent set p(j + A) C C, avoiding (—00, 0]. 
Here z s is taken to be positive for positive real values of z (see |Paz83j page 69, 
|Shu01| pages 83, 94). Positive powers are defined as {~f + A) s = (7 + A) k (~f + A) s - k 
where k £ N and s < k. It can be shown that (7 + A) 3 ^ + A) r = (7 + A) s+r for 
s,r £ M, and that if s £ Z our definition coincides with the usual definition of integer 
powers of (7 + A), see [ShuOlj . In particular we have (7 + A) s = ((7 + A) _s ) _1 . 

We will be interested in powers < a < 1, for which the integral in ([6]) is 
equivalent to 

poo 

{~/ + A)- a = / t a (t + 7 + A)" 1 dt. (7) 

T JO 

as shown in |Paz83j . Section 2.6. This fact will be needed in Section I3~2l when we 
discuss the semigroup theory related to (7 + A) 1 " 13 and (7 + A)' 13 . We will use 
repeatedly in Section [3~5l that the operator (j + A)~ 13 £ C(H S , H S+2 P) is smoothing. 
We denote by £(Hi,H2) the space of bounded linear operators from the Hilbcrt 
space Hx to the Hilbert space H 2 ■ For notational convenience we will sometimes use 
A in place of 1 — j3, in particular when describing the domain H S+2X of (7 + ^4) A = 
(7 + ^) 1 ^ /3 . 

1.4. The Sobolev Gradient. The motivation for equation ([3]) is the desire to 
solve the semilinear elliptic equation 

-Au = V 2 {x lU ) xen (8) 

subject to one of the boundary conditions (0| or (0). Equation ijHJ has an associated 
variational principle. In fact, it is the Euler-Lagrangc equation for the functional 

s (u) = I I (a(x)Vu(x) ■ Vu(x)) + V(x, u(x)) dx. (9) 
Jn 2 

To minimize S, and therefore find a solution to ijSJ), we could consider the steepest 
descent equation 

d t u = -Au-V 2 (x,u). (10) 
The motivation for equation (jTUJ) is that the right-hand side is the unique clement 
g of L 2 so that DS(u)r) = (g,T]) L 2, where DS(u) is the Frechet derivative of S at 
u. This element g £ L 2 is called the L 2 gradient of S with respect to the inner 
product (•, -) L 2. Instead of the standard L 2 inner product, if we used a different 
inner product, we would obtain a different gradient for S. 

We consider the Sobolev space with inner product (u,v)p = ((j + A) /3 u,v)^2 7 
and look for the gradient of S with respect to this space and inner product. That 
is, the unique element g £ such that DS{u)rj = {g,rf)p. We refer to g as the 
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Sobolev gradient of S and write g = V^S^it) (see. |Neu 97 ) . As noted at the end 
of Section 11.21 this gradient depends not only on /3 but also on our choice of inner 
product, which was determined by A and 7. 

We note that in each case, Dirichlct or periodic boundary conditions, we are able 
to use the integration by parts formula 

— / div(a(x)Vu)u dx — / a(x)Vu • Vi> dx, (11) 
Jn Jn 

and we calculate the H ^-gradient as follows: 

DS(u)rj = / a(i)Vw ■ Vr? + V2(x, u)rj dx 



(-div(a(x)Vu) + V 2 (x,u),r]) L 2 = {Au + V 2 (x, u), 77) i2 
((7 + A)' 3 (7 + A)-P(Au + V 2 (x, «)), V ) L2 
((7 + A)~' 3 (7u + Au - 7U + V 2 (x, it)), rj) 
((7 + Af-^u - (7 + A)-f\ lu - V 2 (x, u)), jj)^ . 



Thus, our steepest descent equation in iJ^, 9(U = —X7pS(u), becomes 

d t u = -(7 + Af-^u + (7 + A)-^( 7 u - ^(x, u)), 

which is identical to ((3]). If the solution u(x,t) of ((3|) approaches a critical point, 
that is u(x, t) — > m c (x) as f — > 00, then u c will solve (7+j4) 1_/3 7i c = (j + A)~ l3 ('yu c — 
V 2 (x,u c )), which reduces to |(5J). 

2. Main Result 

We now wish to formulate our main theorem, which is a comparison principle 
for the flow defined by ©. The theorem is actually two theorems, one for each 
type of boundary condition. Thus, in the statement of the theorem, the space 
H s may refer to either of the two types of Sobolev spaces Hq (Dirichlet boundary 
conditions), or Hp (periodic boundary conditions). We will write O to represent the 
space domain, whether it is T d or a smooth, bounded subset of M. d with Dirichlet 
boundary conditions. 

Theorem 2.1. Let V £ C r (fl x R,R), r > 2, and choose 7 such that 7 > 
sup/ a j^gR^xR \V 22 \. Let T > 0, and let u(x,t) and v(x,t) be solutions of 0) for 
t £ [0, T] with initial conditions u(x,0) = uq £ L°° and v(a;,0) = vq £ L°° , where 
the exponent [3 in (0) is taken in the range [3 £ (0, 1). If uq > vq for almost every 
x £ n, then u(x,t) > v(x,t) £ [0,T] and almost every x £ Q. 

Remark. The regularity of V will be the only limit for the regularity of the solution 
u. We show in Proposition 13.61 that u(t, ■) £ H r+S . for any 6 < 2. In particular, if 
r > d/2 — 1 then for u £ L°° we have that the solution u(t, •) £ C°(fi). 

The proof of existence will show that in our case, we have T = 00, which allows 
this semi-flow to be used for finding critical points of S. 

We use the notation L := -(7 + A) 1 ' 13 and X(u) := (7 + A)~ p (ju - V 2 (x, u)) 
so that we can rewrite (|3j) as 

d t u = Lu + X(u), (12) 
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and easily refer to the linear and nonlinear components of the equation as L and 
X. We now briefly outline the strategy for the proof of Theorem 12. II 

We will show that L generates a semigroup and this semigroup satisfies a com- 
parison principle. The theory of semigroups will also allow us to show that, for 
large enough 7, the nonlinear operator X will also satisfy a comparison principle. 
We will then show that solutions to equation (fl"2|) exist for all time, and can be 
expressed via the integral formula 

u(x,t) = e tL u (x) + [ e it - s)L X(u(s,x))ds, (13) 
Jo 

commonly referred to as Duhamel's formula (see |Tay97| , page 272). This formula 
and the comparison principles for e tL and X will allow us to show that the solutions 
to (|12p must also satisfy a comparison principle. 



3. Preliminaries 

3.1. Moser Estimates. The composition V(x,u) will be controlled by the use of 
Moscr estimates for the composition of functions in H s , s G N. If / G C s (M. d x M) 
and <p G H s (M. d ) n L°°(R d ) then 

||/(^)|| s <C s |/|ca(l+|MU). (14) 

We also have that if / G C s+1 and if (f>, ip G H s arc bounded with s G N, then 

\\f(x,</>) - f(x,lP)\\ s < C.|/| . + i(l + |M| S + M\s)H- Tp\\s- (15) 

The constant c s depends on the supremum of cj> and the diameter of f£, see |Mos88j . 
[Mos66| . When proving the existence of solutions to equation ([3]) we will first show 
that they exist in L°° for all time t > 0, and then that they are in fact continuous 
in the domain Q and diffcrentiablc in time. 



3.2. Semigroups and Fractional Powers. We will show that the operator L = 
— (7 + A) 1-13 generates a semigroup and that this semigroup has many nice prop- 
erties, including a comparison principle. Most of these facts will be derived from 
properties of the semigroup generated by —(7 + A). To prove that —(7+ A) does, in 
fact, generate a semigroup, and to begin the investigation of some of its properties, 
we start with the following 

Definition. If H is a Hilbert space and D C H is a dense linear subspace of H, 
then a linear operator B : D — > H is called m-accretive if 

(-Bu,u) > Vu G D, and 

{-B + I)D = H (16) 

The Lumcr-Phillips theorem (see |Paz83] ) states that if B is m-accretive, then 
it generates a strongly continuous semigroup of contractions, e tB . That is, e tB G 
C([0,oo),H) n (^((O.oo),!?), there exists c > such that \\e tB \\ c{H} < e~ ct , and 
if uq G H, then u(t,x) := e tB uo(x) satisfies 

du 

dt = BU 
u(0,x) = uq(x). 
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It is not hard to see that if 7 > then 7 + A is m-accretive on the Hilbcrt space 
H s (with either periodic or Dirichlet boundary conditions) . Let u be an element of 
the dense subspace H s+2 c H s , then 

((7 + A)u, u) s = (( 7 + AY +1 u, u) L * = \\u\\ 2 s+1 > 0. 

Here we have used the inner product on H s as described in Section [L2l 

Standard results from the theory of elliptic boundary value problems also give 
the existence of a solution u £ H s+2 to the equation ((1 + 7)/ + A)u = f, subject 
to either of the two boundary conditions. Thus (1 + 7)/ + A is a surjection from 
H s+2 to H s (this was also discussed in Section ll.2[) . This establishes the second 
condition in (|16p . Hence, by the Lumer-Phillips Theorem, we have established 

Proposition 3.1. For each s > 0, the operator —(7 + A) is m-accretive with 
respect to H s , and therefore generates a semigroup e~ l - 1+A ^ t £ C([0, 00), H s ) (1 
C l ((0,oo),H s ). 

The fractional power (7 + A)~ a for a £ (0, 1) as defined in can be expressed 
in terms of the semigroup e~ ty ~ 1+A ^ t . We have the formula 

1 r°° 

( 1 + A)- a f = — t^e-^+^fdt, (17) 
r W Jo 

as shown in |Paz83j . |Vra03j . The above integral converges in H s for a £ (0,oo) 
because of the bound ||e~ t ^ 7+ ^^|j£(^fa) < e~ Cst , which follows from the fact that 
e -t(7+A) ^ s a sem jg r0U p f contractions on H s . However, (|17l) is equivalent to ([7|) 
only for a £ (0, 1). The representation above will be useful in proving a comparison 
principle for the operator (7 + A)~@ in the next section. 
A theorem in |Kat60j gives 

Proposition 3.2. For each s > 0, the operator L := — (-f + A) 1 "^ is m-accretive in 
H 3 , and therefore generates the semigroup e tL £ C([0, 00), H s ) PI C 1 ((0, 00), H s ). 

Instead of referring to Kato, one can see this fact in the self-adjoint case by an 
argument similar to the one for — (7 + A). In particular, ((1 + 7)/ + A) a maps 
jjs+2a on f fjs^ as discussed in Section [L2l And for u £ H s+a we have 

<( 7 + A) a u, u) s = ((7 + A) s+a u, u) L , = \\uf s+a > 0. 

3.3. Further Estimates for e tL . If H is a Hilbert space and B is self-adjoint and 
m-accretive on H (this implies B is regular m-accretive, as defined on page 22 of 
|Sho97j ). then for any integer n > 1 and any u £ H, we have e tB u £ D{B n ) and 
that 

^^s)"' (18) 

see page 29 of |Sho97j . This result applies to fractional powers of —B. In fact, if 
B is as above, and a £ (0, 1) then there exists a constant C Qi t such that for any 
t £ (0, T] one has 

\\(-Bre tB \\ c{H) <C a , T (j) a . (19) 

Furthermore, if we set Y = D((—B) a ), the domain of (—B) a , endowed with the 
graph norm ||uj|y = \\u\\h + ||(— B) a u\\u , then 

\\e tB \\c(H,Y)<C a , T (jy, and 

\\e tB -I\\c(YM)<C' a ^. 
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For further details on (HI), dHJ), and {20]) see Section 4.1 of |dlLV09j . 

We can use these facts to show that e tL is a smoothing operator and to establish 
operator bounds on e tL and e tL — I. We say that e tL is smoothing if it increases 
the regularity of a function as measured by membership in different spaces. That 
is, an operator is smoothing if it maps a space of functions into another space of 
smoother functions. In our case, if u 6 L 2 , then for any s > 0, e tL € H 3 . To see 
this, we first note that because L is sclf-adjoint and m-accretive on H s for s > 0, 
estimate ([!%]) gives 

<(£)'■ (2D 

Recall, as we established in Section 11.21 that L : H S+2X — > H s and that the inner 
product on H a is given by (u,v) s = ((7 + A) s u,v) L 2 = ((—L) s / x u,v) L 2, where we 
set A = 1 — p. Note that A > 0. If u 6 H s , we compute 

\\e tL u\\ 2 s+2nX = (e tL U,e tL u) s+2nX = ((-Lf nX / x e tL U,e tL ) s 

2n 



((-L) n e tL u, (—L) n e tL ) s = (-l) 2n \\L n e tL u\\ 2 s < 
Thus, e tL is smoothing in the sense described above, and we have the bound 

\\e tL \\c(H*,H>+^) < {jj^J ■ ( 22 ) 

To apply in the case B = L and H = H s , we have that Y = H s+2aX . Then 
estimates (0 and (gHl) yield 

\\{-LTe} L \\ c{HS) <C a , T (^) a , (23) 



l \\c(h;h-+™ ) <C a . T (^y, (24) 

,\e- -I\\ c( H^, Ms) <C'^ T t a . (25) 
In fact, can be obtained from a calculation similar to that above for \\e tL 



,t£i 

— J !l£(ff s + 2 ° A ,ff=) ^ L/ Q! ,T'" • 

3m a calculation similar to that above for \\e tL u\\ s+2nX - 
In particular 

j) HI- 

Thus, for a € (0,1), He^lU^.^-A) < C Q , T (?)". 

3.4. Subordination and Comparison. The operator 7 + ^4 is uniformly elliptic 
so the maximum principle for parabolic equations applies to 
du 

— = -{ 1 + A)u, u(0,x)=u (x). (26) 

Thus, a solution to (|26[) on the interval [0, T] obtains is maximum on the boundary 
of n x (0, T], but not at fi x {T}, (see [PW67]). Therefore, if u > 0, then for t > 
we have e - *^ 7 "*"" 4 ^ > 0. We will use this fact to establish comparison principles for 
X and e tL . We begin with 

Proposition 3.3. Let 7 > sup^ y {V22(x, y)}. Then X satisfies a comparison 
principle. That is, if u > v a.e., then X(u) = (7 + A)~^(-fU - V 2 (x,u)) > 
{j + A)-f 3 (ju-V 2 (x,u))=X(v) a.e. 
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Proof. Since 7 > sap xy {V22(x,y)} ) then ju — V2(x,u) is increasing in u. Thus, 
u > v implies 71* — V%(x, u) > "fv — V%(x, v). Then, as mentioned above, the 
maximum principle for parabolic equations implies 

e - t ^+ A \ 1 u-V 2 {x,u)) > e- t ^ +A) { 1 v-V 2 {x,v)). 

Hence, for t > 0, 

thus the representation of (7 + A)~@ in equation (fTT|) implies that (7 + A)~@("fu — 
V 2 (x, u)) > (7 + A)~P(^v - V^x, v)), and we conclude that X(u) > X(v). □ 

To show a comparison principle for e tL , we again use the semigroup generated 
by 7 + A. However, we will need a different representation of e tL , which is given 
by the subordination identity of Bochncr |Boc49j . For a > 0, t > 0, r > 0, and 
< a < 1 we define 

1 />cr+ioo 

t ,a(r) = — / e T2 ~ tz °dz 

27T* ,/o-ioo 

and t , Q (-r) = if r < 0. Additionally, <j> tt0l (r) > for all r > 0. Then, for B 
m-accretive, we can represent e~^~ B ^ as 

/•oo 

e-* ( - sr =/ e TB t . Q (r)dr, t > 0, (27) 
Jo 

as is shown in |Yos74j Section IX. 11. This provides a representation of e tL as 

/•OO 

e tL = e -tfr+A)* = e -r( 7 +A) 0tA(r)dTi yt>0j (28) 

Jo 

with (f>t,\(r) > 0, which will be key tool in establishing 

Proposition 3.4. If u > v a.e. in il, then e tL u > e tL v in fi. 

Proof. If it > u a.e. then for each t > and r > 0, we have 

e-^ +A ^ tA (r)u > er^ +A ^ t . x {T)v 

because 4>t.\ > and e~*' 7+j4 ' satisfies a comparison principle as explained above. 
Integrating both sides of the inequality yields e tL u > e tL v. The fact that e tL 
is smoothing, in the sense of inequality (|2"4"j) . guarantees that e tL u and e tL v are 
continuous from fl to M, and therefore e tL u > e tL v for all x £ O. □ 

In preparation for the proof of existence and uniqueness of solutions to (|12D . 
which will require L°° estimates on X and e tL , we will show that X and e* L are, in 
fact, locally bounded maps from L°° to itself. This is clear for e tL by the remark 
at the end of Section [3721 because for t > 0, e tL u <G H s for arbitrarily large s, and 
therefore it is in L°° by the Sobolev embedding theorem. However, we can use the 
comparison principles for X and e tL to provide explicit bounds. 

Proposition 3.5. X : L°° — > L°° is locally bounded with \\X(u)\\l°° < 7 A |I u I|l oc + 
7 || V2IIL00 • Additionally, for each t > 0, e tL : L°° — > L°° is a bounded linear map 
and\\e tL \\ c{L ^<e-~< Xt . 
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Proof. Let u G L°° , and set u — ||w|ji=o, so that u < u. Then the comparison 
principles established in Propositions 13.31 and 13.41 imply that X(u) < Xiu) and 
e tL u < e tL u. Hence, if X and e tL are bounded on constant functions then they are 
bounded on L°° . 

Consider first X, and set Cy — yu + ||V2||z,°°. We see that the boundedness 
of V% gives yu, — V^ix^u) < Cy. The integral representation of (7 + A)~@ in 
equation (| 1T[) and the same reasoning as in the proof for Proposition 13.31 imply 
X(u) = (7 + AyP^u - V 2 (x, u)) < (7 + AY^Cy. Thus, to establish bounds on 
X , we only need to understand how (7 + A)~@ acts on constants. 

A consequence of Proposition 10.3 from [ShuOlj . page 93, is that if ip is an 
eigenfunction of 7+ A with eigenvalue /i, then ip is also an eigenfunction of (-f+A)^ 13 
with eigenvalue . But (7+ A)Cy = jCy because A is a second-order differential 
operator. Hence (7 + A)~^Cy = 7 _/3 Cy and 

X(u) < X(u) < 7 - /3 ( 7 S+ \\V 2 \\l°°) = 7 A |M|l~ +J~ I3 \\V2\\l*, 

establishing the first claim in Proposition 13.51 

To bound e tL we examine how L acts on constants. Using (7 + A)~^u = 7 _/3 u, 
we calculate 

Lu = -(7 + A) x u = -(7 + A)(j + Ay fi u = (7 + Ay~ p u = -j^^u = -~j x u. 

Now e tL is the semigroup generated by L, so we know e tL u solves dtu = Lu, 
u(0) = u. But Le tL = e tL L, hence 

d t e tL u = e tL Lu = e tL (-j x u) = -y X e tL u. 

Thus e tL u solves dtu = — j x u, u(0) = u. Hence e tL u = e |,A *u, and we have 
e tL u < e tL u = e'^u. □ 

3.5. Existence and Uniqueness of Solutions to Equation (112[) . The com- 
parison result in Theorem 12.11 requires only the existence of the flow generated by 
(TT2"]) for some short time T > 0. However, the motivation for studying this flow is 
to find critical points of the functional @ , for which the flow must be defined on 
all of (0, 00). In this section we establish the following 

Proposition 3.6. // the potential V G C r+1 (fi x M,R), r > 1, and u Q G L°°(n), 
then for every S < 2 there exists a unique solution, u(x,t) G C((0, 00), H r+s r\L°°), 
to {23> with u(x,0) = u (x). Ifr > 2 then u(x,t) G C7 1 ((0, 00), H r ~ 2X ). 

We begin by showing existence and uniqueness of a mild solution in L°°, and 
use the smoothing properties of the flow to obtain the desired regularity. 

Definition. We say u is a mild solution of the equation dtu = Lu + X(u), with 
w(0, x) = uq(x) if u satisfies 

u{t,x) = e tL u {x) + [ e^- T)L X{u{T,x))dT. 
Jo 

We consider the map 4", on C([0, T], L°°) : u(0, x) = u G L°°} given by 

Vu(t) = e tL u + f e(*- r)i X(u(T)) dr. 
Jo 

The mild solution, u(x, t), of (Tj"2"|) will be the fixed point of 'J. 
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Lemma 3.7. For T g R, define W T = {u £ C([0, T], : u(0,ar) = w g 

L°°}. T/ien for small enough T > 0, \& is a contraction on Wt- The size of T is 
independent of Uq. 

Proof. The norm on W will be denned as ||u||oo,T = su P r e[o,T] ll u ( T )lk 00 - For any 
t € [0, T], it is clear that ||w|joc.t < ||w||oo,r- Let u g Wt, then \I>u(0) = uq, and we 
have 

||*«(*)||l- < ||uo||l-+ f e-^ {t - T) (l X \\u(T)\\ L ~ +1 -P\\V 2 (x,u(x,T))\\ L ~)dT 

Jo 

< HuolU- + (1 - e-^ t )(7 A ||w||oo,T + T^II^IU-) < Or, 

where CV is finite for finite T. This follows directly from Proposition 13.51 

Note that by the differentiability assumptions on V, we know that for any x £ f2 
and any y x ,y 2 £ R, |V 2 (x,?/i) - V 2 (x,y 2 )\ < \V 2 \ c i\yi ~ 2/2 1 - To see that * is a 
contraction, we compute 

\\y u (t)-yv(t)\\ L ~ < f \\e^ L (X(u)-X(v))\\ L ^dr 

Jo 

< f e-^( t -^\\ 1 x {u~v)+ 1 -P{V 2 {x,v)-V 2 {x,u))\\ L ~dT 
Jo 

< (1 - e-* 1 )^ + T^l^aloOHu - "IUt = Ct||u - vIIoo.t, 
where C = 7 A (7 A + 7~^|V r 2|c' 1 ) depends only on 7, j3, and V. Here we have used 
only for convenience the fact that 1 — e -7 * < 7 A t for t > 0. Choosing T = 
ensures that Ij'J'M— ^^Hoo.t < w|| 00, t, and thus ^ is a contraction on Wr- This 
choice of T depends depends only on 7, /3, and V. In particular. T is independent 
of the initial condition uq £ L°° . □ 

Proof of Provosition \3.6[ Lemma [3~7l ensures that, for T sufficiently small, ^ has a 
unique fixed point u £ C([0,T],L°°) satisfying 

u(t) = e tL u + f e^- T)L X(u(T))dT, 
Jo 

establishing the existence and uniqueness of a mild solution to (fl2|) in L°° . Because 
T was chosen independently of uq, we will also have existence on [0, T] with initial 
condition u(x,T), which gives existence on the interval [0, 2T] for initial condition 
uq. This can be repeated indefinitely and we have L°° existence on [0,oo). 

The operator (7 + A)^ 13 is smoothing in the sense that it maps H s into H S+2 P. 
Specifically, for w £ H s , 

\\{l + A)-> 3 w\\l +w = (( 1 + A)-P W ,( 7 + A)-ew) s+20 

= ((7 + A)- 2p w, w) s+2p = (w, w) s = \\w\\ 2 s , 

which gives \\(-f + A)~P\\ C ( H s_ HS +2f3- ) = 1. This implies that X will have a smoothing 
property too. However, this will be limited by the regularity of V. Recall that f2 
is bounded, so L°° cL 2 = H°, hence (7 + A)~P maps L°° into H 213 , and for any 
u £ L°°, and 

\\X(u)\\ 2 p<\\ 1 u+V 2 (x,u)\\ < MMl* + \\V 2 (x,u)\\ L 2 < |fi| 1/3 (7Nl£~+l|V 2 ||L~). 

Thus X(u(t)) is bounded in H 2 " as long as u{t) is bounded in L°°. Note that for 
any fixed t £ [0, 00) we know that ^(r)^^ is bounded, where the norm || • \\oo,t 
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is the same notation as in the proof of Lemma (|3.7|) . To see that u(t) is actually a 
solution in H 2 " , we simply compute, for a fixed t > 0, 

IK*)I|20 < 

< 
< 

From here we wish to repeat this process to show that u(t) is actually a solution 
in H At3 . However, this will require a bound on the composition V(x, u(x, t)) in the 
space H 2 @ . For this we will need to employ the Moser estimates (HU), but these 
estimates only apply for H n , n G N. This difficulty can be handled by splitting our 
analysis into two cases, first for (3 G [1/2, 1) and then for (3 G (0, 1/2). Recall that 
e tL is smoothing in the sense of estimate (|22|) . so that if uo G L°° C L 2 then for all 
t > 0, e tL ua G H s for any s > 0. Thus, the term ||e* L ito||n+i that appears in the 
estimates below will not impede to our regularity-building scheme. 

Suppose (3 G [1/2,1), then the solution u(t) G L°° of {12]) is bounded in H 2f) , 
but 2/3 > 1, so < \\u\\ 2f} . Thus, u(t) is an H 1 -solution of JT5J). To show that 

ii G iJ", we use induction on n, with n = 1 just established. Now assume w G if" 
with n < r, then by JTJ| we have ||V 2 (x,M)|| n < cy{\ + ||u|| n ), where cy depends 
on IMIoo.t and V. Because 2/3 > 1 we have, for any w G H n , + A)~^u;||„ +1 < 
HHU+i-2/3 < \\ w \\n- Using these facts, wc compute 

r t 

||«(*)IU+i < 
< 

< 

< 
< 

re[0,t] 

Thus, u{t) G H n+1 , for any n < r. Improved regularity to u(t) G £T r+ * for any 
5 < 2 will be shown after the case /3 < 1/2. 

Between lines one and two of the above computation we have used the fact that 
e tL is a contraction semigroup on H n+1 , i.e. \\e( t ~ T * >L \\c(H™+ 1 ) < 1- We did not, 
however, take advantage of the smoothing property of e tL . This is something we 
cannot afford to waste when f3 G (0, 1/2). 

Suppose 13 G (0,1/2), then A = 1 - /3 G (1/2,1). Let u(t) be the L°°-solution 
of (fT2")l . The previous method of using the smoothing properties of (7 + A)^ 13 
fails because the smoothing factor of 2/3 is too small. Fortunately, this is precisely 
when estimate provides a large smoothing factor from e tL . Applying (f!M)) with 
a = 1/2A gives 

\\e tL \\c(H»,H^) <C x , T t- 1/2X . 



\\e tL u \\ w + / \\e^ L X{u{r))\\ 2p dT 
Jo 

||e"^o|| 2 /3 + M 1/2 / e-^(*-) (7||«(r)|Uoc + ||y a |U-)dT 
Jo 

\\e tL uo\\ 2 + |^| 1/2 ^ (l - (7||tt||ao,t + ||^|U~) < 00. 



e" J W ||n+l+/ ||(7 + ^) IJ {JU{T) -V 2 (x,u(T)))\\ n+1 dT 

Jo 

e tL u Q \\ n+1 + [ \\~/u(t) - V 2 (x,u(T))\\ n dT 
Jo 

e tL u Q \\ n+ i + / 7||u(T)|| n + (v(H-||«(T)||„)dT 
Jo 

e tL u |L+i + Ct sup (1 + ||u(t)|L) < 00. 
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The cost of the smoothing is the factor of tr 1 / 2X . However, 2A > 1 so this is 
integrable on [0,T], so we can proceed with a similar argument to the case for 
(3 G [1/2, 1). Just as we did not use the smoothing properties of e tL in the previous 
case, we do not need the smoothing properties of (7 + A)^ 13 in this case. To begin 
induction on n, we note that u(t) G L°° C H°. If u(t) G H n with n < r, then 

IK*)l|n+l < 
< 
< 

where C depends on cy, 7, A, and T. Thus, u(t) G H n for all n < r + 1. 

In either of the two cases, j3 > 1/2 and (3 < 1/2, we have established that u G H r . 
From here we can easily improve to u G H r+2 ^ because 

\\u{t)\\ r +2i3 < ||e tL u ||r+2/3+ / l\\u{r) \\ r + cy(l + \\u(t) \\ r )dr, 

Jo 

regardless of the size of (3. Now we can push a little further and consider \\u(t) || r _|_2/3+e 
for some e > 0. We can use estimate (fiM)) to obtain 

\\e^ L X( U (T))\\ H r + ^ a <C*(t- r)-«(l + |Kr)|| r ), 

provided a G (0,1). The right side of the inequality is integrable on [0,£]. Thus 
u(t) G H r+2fj+2Xa for any a < 1, which gives u(t) G H r+S for any 5 < 2(3 + 2A = 
2/3 + 2(1 - (3) = 2. Hence u G C([0, T],H r+s n L°°) for any 8 < 2. 

To finish the proof of Proposition ^. 61 we must show that the derivative, u tl of u 
exists in H r ~ 2X so that u is actually C 1 and not just a mild solution. To see this, 
set p = r — 2A, and consider Rh = \ \\u{t + h) — u(t) — hu t \\ p , where 

u t :=Le tL u + f Le (t ~ T)L X(u{t)) dr + X(u(t)), 
Jo 

Showing Rh — > as h — > will prove the desired result. We have 
Rh < \\^(e hL ~I)e tL u -Le tL u \\ p + 

+ II f ({(e hL -I)-L) e^ L X(u(r)) dr\\ p + 
+ II \ e^ h ^ L X{u{r)) dr X(u(t))\\ p 

For brevity, we will refer to the three terms on the right-hand side of the inequality 
as h,l2, and I3. I\ goes to zero with h because L is the generator of the Co- 
semigroup e tL on H p and the fact that e tL Uo G D(L) for all t > 0. In I2, we can 



\\e tL u \\ n+1 + / \\e^ L X(u(r))\\ n+1 dT 
Jo 

||e ti u ||n+i+ / C x .Tt- 1 / 2X \\( 1 u(T)~V 2 (x,u(T))\\ n dr 



\\e tL U \\n+l+C SUp (l+||«(T)||„)t 1 - 1 / 2A <C50 J 
re[0,t] 
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replace (e hL — I) by J Q Le aL da and we see that 



1 



h = ||/ e^~ T > L - / Le" L da-L)X(u(r))dT\\ p 



< 



I e(t ~ T)L T I {e aL - I) LX(u(T))dadT\\ p 
Jo Jo 

e ° L - I \\c(H'P+ 2 '-^,HP)\\ LX ( u ( T ))\\p+'2/3\dadT 



t 2 rh 



a h J 



< I \ I C 2 ^da\\LX(u(r))\\ 



lo h Jo 

< C' 2 h [\\\(u(T))\\ r +Cv(\\u(r)\\ 00 )(l + \\u(r)\\ r )) dr. 



We have used \\e^ t ~ T ' )L \\c(Hp) < 1 and, between lines three and four, applied ([25]) 
with a = (3. The estimate in the final line of the calculation follows from the bound 
||(7 + ^4) 1_2/3 w|| s < ||HU+2(i-2/3), which applies to the operator LX(u). Hence, 
the fact that p + 2(3X + 2(1 - 2(3) = r - 2/3 2 < r allows the use of (M]). We know 
that for any fixed t > 0, ||u(t)||,. and HutY)^ are bounded on t £ [0,t], so I 2 can 
be made arbitrarily small for a suitable choice of h. 

Without loss of generality, we assume that h < 1 and we set Mi = max t < T < t+1 ||u(t)|| 
and M 2 = max t < r < t+ i ||u(t)|| p . Let e > 0, by continuity there exists h G (0, 1) such 
that \\u(t) — u(t)\\ p < e for |r — t\ < h. Thus we have 

1 rt+h 

I 3 < - \\e^ h -^ L (X(u(r))-X(u(t)))\\ p dr 



1 



t+h 



< I I h{u{r)-u{t))-(V 2 {x,u{T))-V 2 {x,u{t)))\\ p ^dT 



1 



t+h 



< -J C v (Mi)(l + 2M 2 )\\u(r) - u(t)|| r _ a 

< ^f- f + \Ht) - u(*)|| r _a < C' VM e l -h < eC' VtM , 

so that ^3 can be made arbitrarily small. We have again used ||e^ + ' 1 ^ T ' i j|£(^ P ) < 1 
and that p — 2(3 = r — 2. We have also used (fTS)) between lines two and three above. 
Thus we have shown I\, I 2l I3 — > as h — > 0, establishing u £ C 1 ([0, 00), H r ~ 2X ), 
completing the proof of Proposition 13. 61 □ 



4. Proof of Theorem 12.11 

We have established comparison principles for the semigroup e tL and the op- 
erator X, as well as the existence of solutions to (fl"2)) . To emphasize the initial 
conditions, it will be convenient to write the solutions of (|12p as u(x, t) = &tUo- 
Hence, we aim to show that if «o, v o € L°° an d uq > Vq then $ t uo > &tVo on a 
short time interval [0, T]. This will follow from the iteration method below. 
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4.1. Iteration Method. For u G L°° we define F^u = e tL u, and the jth iterate 
of u as 



F t J+1 u = e ln u + 



Jt-r)L 



X(F^u)dT, 



denned on some interval [0, T]. Fj! + is well defined because X and e tL are both 
bounded maps from L°° to itself, as shown in Proposition 13.51 



Proposition 4.1. LetT>0. If u > v then F ] t u > F{v for all t G [0,T]. 

Proof. Assume u > v. Then by Proposition 13.41 F t °u = e tL u > e tL v — F®v. We 
assume that F^u > -F/y and proceed by induction on j. By Proposition 13 -31 we 
have that X(F^u) — X(Ffv) > 0. Once again invoking Proposition 13.41 we have 
e (t-s)L[X(Fj>u) - X(Fiv)] > 0. Hence 



Fl +1 u-F^ L v = e th {u- 



)) + f e^- 3 ^ [X(Fj>u) - X(F*v)] ds>0 
Jo 

because the integrand is positive and e tL (u — v) > 0. 



□ 



We need to show that this iteration converges to the solution in Theorem l2.11 so 
we now focus on a single initial condition. u . For notational convenience we write 
w J '(t) in place of F^uq, and u°(t) = e tL UQ. Thus for each j G N we have 



U i +1 (t) =e tL u + / e (t - s)L X(u^s))ds. 
Jo 

Proposition 4.2. If uq G L°° then there exists aT > such that for all t G (0, T], 

u j (t) G L°° n -ff s , /or a// < s < r + 1 and every j G N. 



Proof. A consequence of Proposition 13.51 is that e tL and X are bounded on L c 
We use this to compute 



\W +1 (t)\\ 



< 
< 
< 
< 
< 



e tL u \\oc + 



"olloo + 7 1 max H^Hoo 

0<.s-<T 

^o| oo 

+ T(|| 

^0|| OO 

+ T max Hoc) < . . . 

0<s<T 

wo||ooT 1 + HuollooT 2 + . . . + \\u \\ooT j + T3 +1 \\u Q \\ 



1 _ T'J+ l 
l-T 



+ T3+ 1 \\u \\ oc , < \\u \ 



l-T 



= C 



Where, without loss of generality, we have assumed T < 1. So we know that 
each iterate u? is contained in the ball with radius Co in L°° for any j and for all 
t G [0, T]. The i°° bounds allow the use of the Moser estimates (fT4")) . For integer 
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k < r + 1 we have 

\\u j+1 (t)\\ k < \\e tL u \\ k + [ \\e^ L X( U i(s))\\ k ds 

Jo 

< \\e tL u \\k + C T max (1 + \\u\s)\\ k ) 

0<s<T 

< \\e tL u \\ k + C T(l + \\u Q \\ k + C T max (1 + ||u , '- 1 («)||k)) < ... 

0<s<T 

< ||e tL Uo ||fe + C T(1 + \\u \\ k ) + ... + {C Ty(l + \\u \\ k ) + (C T) j+1 (l + \\u Q \\ k ) 

< sup \\e^u \\ k+ { \ + ^l k) =C k . 

0<t<T -L — OoJ 

□ 

We have assumed that TCq < 1 and T < 1 . which require only that T < 1+ ||„ || — • 

Proposition 4.3. u J converges in C([0, T], H r ). 
Proof. 

\\u j+1 (t) - u\t)\\ r < f \\e^ L [X(u j (s)) - Xiu^isMlr ds 
Jo 

< t max \\X(u j (s)) - Xfu^is^lL 

0<s<t 

< t max C\V\ r C (l + \\u j (s)\\ r + !|^- 1 ( S )|| r )||^'( S ) - u^ 1 ^ 

0<s<t 

< t max C*\\u j (s) - n^VsllL, 

~~ 0<s<t 

where C* depends on \V\ r , Co, and C r . So we have 

max \\u j+1 (t) - u j (t)\\ r < C*T max ||u J (s) - u^^sML. 

0<t<T 0<t<T 

We can set T = -^F, which will ensure that the sequence w-'(i) is Cauchy in 
H r . The argument for this is the following. First, notice that if — u 3 \\ < 

- ui-i-W, then \\u j+1 - u j \\ < {\y\\u l - it°||. So for any e > 0, choose N such 
that — u°\\ < e. Then we have, for m > n> N, 

\\u m - u n \\ < -\\u m - u m - l \\ + -\\u' n - 1 - u m - 2 \\ + ... + -\\u n+1 - u n \\ 



< 



2" 11 2" 2' 



\ 7TC— 1 / 1 \ m ~ ^ ( \ 



2/ \2J \2 

m — 1 



■>m— |j „, 1 „,0 1 



< ( - J [l + 2 + 2 2 + ... + 2 r " " -j || 
1^ 2 m - n \\u x -u a \\<(l\ 



□ 



Corollary 4.4. There exists aT > such that $ t satisfies a comparison principle 
on the interval [0,T]. 
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Proof. From Proposition 14.31 we have the existence of u°° <E C([0, T], H r ) such that 
v? — > u°° in C([0, T],H r ). This function u°°(f, x) = lim^oo F^uq(x) must satisfy 

u°°(t,a;) =e* L w + / e ( *- s)i X(ii°°( S , a;)) da, 
Jo 

and therefore, by Proposition 13. 6i u°°(t,x) = &tUo- By Proposition 14. 11 if uq > vq 
then F^uq > F^vq and therefore $(Uo > $(Uo on [0, T]. Thus we know <I> t obeys a 
comparison principle on a small time interval [0, T]. □ 

This establishes the comparison principle on a finite time interval [0, T] and 
therefore concludes the proof of Theorem 12.11 To see that this comparison holds 
for all time t > 0, we have the following lemma. 

Lemma 4.5. If uq > vq for a.e. x G Q,, and there exist a time t\ > for 
which <& tl uo(x) < ^avo^x) on a set of positive measure then for every t > 0, 
<fr t Uo( x ) < < l 5 tWo(a;) on a set of positive measure. 

Proof. Let Ua(x) > ^0(2^) for a.e. a; G O. Suppose there is a first time t\ such that 
We define the family of $tUo(a;) < &tVo(x) on a set of positive measure. Then on 
this set of positive measure 

e tlL (u (x) -v (x))+ / e^- T > L (X($ T u (x)) - X(<S> T v Q (x))) dr < 0. (29) 
Jo 

However, e tlL (uo(x) — vq(x)) > for all x G O by Proposition 13.41 We have by 
assumption that for all r G [0, £i) $ r w(x) > $ r u(x) for a.e. a; G O. Hence, by 
Proposition DL3] for all r G [0,<i) X($ T u(y)) - X{<f> T v(y)) > for a.e. x £ fl 
Again applying Proposition 13.41 we have for any, x € 

X{$ T v{y))) > for all r € [0,ti). Thus for each r G [0,ti) the integrand in j29]) 
is a non-negative function on f2, and therefore the integral (which is a function in 
H r+1 ) must be non-negative for a.e. x G f2. Therefore the left side of inequality 
(|29[) is the sum of the two terms that are non-negative for a.e. x and cannot be 
strictly negative on a set of positive measure. □ 

Combining Corollary 14.41 and Lemma l4~5l we have that <J> 4 satisfies a comparison 
principle on the interval [0, T], for T > and therefore $ t satisfies a comparison 
principle on [0, 00). 

5. Fractional Elliptic Equations 

We can extend the methods above to gradient descent equations for energy func- 
tionals of the form 

S a (u) = (u, A a u) L 2 + / V{x,u)dx, 
Jn 

where A is given by ^ and a G (0, 1). The Euler-Lagrange equation for S a is 
A a u + V 2 {x,u) = 0. 

For a fixed a > 0, we use the inner product for the Sobolev space H a given by 

(«,«>q = {{l + A a )u,v) L 2 
we have the inner product on H ar given by 

(u,v)* as = ((7 + A a yu,v) L *. 
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Now consider a £ (0,1). To calculate the Sobolev gradient of S a , we first note 
that the derivative of S a is DS a (u)r] = (rj,A a u + V 2 (x,u)) L 2. Thus, the Sobolev 
gradient of S a in H a/3 , (i £ (0, 1), is calculated as 

DS a (u) V = (v,A a u + V 2 (x,u)) L 2 

= (r), ( 7 + A a )P(n + A a )~ l3 (A a u + 1U - 1U + V 2 (x, u))) L 2 

= (77, ( 7 + A a yP(A a u + 1 u- lu + V 2 (x, u)))* ap 

= ( V , (7 + A a ) 1 ~ l3 u - ( 7 + A a )-P{ lu - V 2 (x, uW af3 . 

Hence, the gradient descent equation is 

d t u = -(7 + A a f- fi u + (7 + A a )- f \ lu - V 2 (x, it)). (30) 

More concisely, we write dtu = Lu + X(u) with Lu := — (7 + A a ) 1 ~ l3 u and X(u) := 
{ 1 + A a )-^(ju-V 2 (x,u)). 

Recall that the maximum principle for parabolic equations ensures that the semi- 
group e~ tA satisfies a comparison principle. The Bochner subordination identity 
(|27|) allows us to write the semigroup generated by —A a as 

/•oo 

e _M ° = / e- rA <f>t, a { T ) dT i *>0. (31) 
Jo 

Thus, by the same argument as in the proof of Proposition 13.41 the comparison 
principle for e tA guarantees that e~ tA satisfies a comparison principle as well. 

Analogously, — (j + A a ) generates the semigroup e~ l ^ 1+A \ which also satisfies 
a comparison principle. Just as in (|17|) we can define the real powers as 

r(p) Jo 

and positive powers are again the inverses of negative powers. Thus, because we 
have a comparison principle for e~*' 7+A °' we have that (7 + A a )~P also satisfies a 
comparison principle. An argument as in the proof of Proposition 13.31 shows that 
if 7 > sup.,, ^ |V22(x, then X will satisfy a comparison principle. 
Employing the subordination identity once again, we have 

poo 

Jo 

Hence, we see e tL satisfies a comparison principle because 4>t,i-p{T) > for all t > 
and e~ T ^ +A ' satisfies comparison principle, just as in the proof of Proposition l3.4l 

Smoothing estimates like (|24|) and (|25|) follow for e tL just as they did for e tL 
in Section [321 The existence of solutions for dtu = Lu + X(u) follows from these 
smoothing estimates and the arguments from Section [3.51 

Finally, the proof for a comparison principle for the flow defined by dtu = Lu + 
X(u) follows from the Duhamcl formula 

u(t,x) = e~ tL u a (x) + [ e- {t ~ T)L X{u{t, x)) dr, 



and the iteration argument from Section [U 
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5.1. Constant Coefficients. In the case of periodic boundary conditions and if 
the matrix a(x) is constant (i.e. independent of x), the situation is simplified 
because we can write down concrete formulae for the operators, and we can use 
classical Fourier analysis in place of some of the abstract semigroup theory. For 
instance, we could have avoided the theory of fractional powers of operators and 
instead used the definition (7 + A) a u = ((7 + 47r 2 £ T a£) Q u(£)) v 

One can also write down the semigroup e~ t ^ +A ^ > = e -t ( 7 ~ A ) = e~ 7 *e tA . The 
operator e tA is convolution with the heat kernel, written as 

(*) - jAln f E e- lx - y+kl2/4 Hy) dy. (32) 

Combining this with equation (|28[) yields the comparison principle for e tL immedi- 
ately. 

6. Application to Aubry-Mather Theory for PDEs 

We can use the comparison principle from Theorem 12.11 and the results from 
jdlLV09| to further develop the Aubry-Mather theory for PDEs to the case of a 
general elliptic problem of the form (JSj). For this section, we restrict discussion to 
periodic boundary conditions. We require the potential function V and the matrix 
coefficient functions a^(x) to be periodic over the integers. That is, 

V(x + e,y + l) = V(x,y) V(e, Z) e Z d x Z, V(x, y) G R d x E 

a lJ (x + e) = a 11 (x) Ve G Z d , Vx G R d 

which we write as V : T d x T — > R and a u : T d — > R. An important class of 
functions is given by the following 

Definition. A function u : R d -> K is said to have the Birkhoff property (or u is a 
Birkhoff function) if for any fixed e € Z d , and I G Z, u(a; + e) — (m(.t) + Z) does not 
change sign with x. That is, u(x + e) — (u(x) + I) is cither > or < depending 
on the choices of e and I, but not x. Any such function can be seen as a surface in 
T d+1 without any self-crossings. 

For a fixed oj € R d , we define B u = {u : R d -> R|u is Birkoff, u{x) — U) ■ x G 
L°°(IR <i )}. We note that u G S w if and only if u is a Birkhoff function and if for 
any e G Z d and Z G Z, cither it(x + e) — u(x) — I < or > according to whether 
w-e — Z<0or>0 (see [Mos86j ) . The vector oj is referred to as the frequency 
or the rotation vector of the function u, and is a natural generalization of the one 
dimensional notion of rotation number. We have the following result. 

Theorem 6.1. Let V £ C 2 (T d x T,R), and let A be a self-adjoint, uniformly 
elliptic operator given by Au = -div(a(x)Vu) with coefficients G C°°(T d ,]R). 
Then for any uj G M. d , there exists a solution u G B u to equation (0). That is 
Au + V 2 {x, u) = 0, u{x) - LJ-xG L°°(R d ), and u is Birkhoff. 

The method of proof from !dlLV09J is to first show the result holds for rational 
frequencies (i.e. for any ujn G jjZ d with JV e N). Then we obtain solutions for 
arbitrary oj G M. d as limits of solutions with rational frequencies. Passing to the 
limit requires an oscillation lemma of De Giorgi-Moser type, as well as classical C e 
elliptic estimates. 
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For minimal solutions, the oscillation lemma was given by Moser in Theorem 
2.2 of |Mos 86j. This gives a bound on the supremum of the solution un by the 
norm of its associated frequency lon <e jj^ d , and independent of N. We provide 
the arguments for the oscillation lemma in Section [6.21 

To have u>n — ► w, we must have N — > oo, and the size of the fundamental 
domain. NT d , becomes unbounded. If the bound on un depended on the size of 
the domain, then the C e -estimates on Vun would degenerate. Instead, because 
the bounds on the un are in terms of lun, the C e -estimates on Vmat are uniform in 
N. Hence we can conclude the convergence of iijy to a continuous function u with 
associated frequency ui. 

Note that the regularity assumption on V is weaker than in Theorcm l2.1l bccause 
we can choose a smooth initial condition (i.e. u (x) = to ■ x). 

6.1. Rational Frequencies. 

Lemma 6.2. If w G -^Z d and u(x,i) solves f3|) on NT d with initial condition 
uo(x) = lo ■ x, then there exists a sequence t n — > oo as n — > oo such that u(x, t n ) — > 
it* in L 2 and u* solves 

Proof. Suppose u(x,t) solves ([3]) on [0,7V] d with periodic boundary conditions and 
intial condition u Q = ui-x. In Scction[TT2]we showed that DSn(u)i] = ((-f+A) 1 ~^u— 
{^+A)~ 13 (p/u — V2(:r, u)) , rj)p. This together with the assumption that u(x, t) solves 
(J3j) allows us to show that Sjv(w(x, t)) is decreasing in t. More precisely, 

^-S N (u(t)) = DS N (u(t))d t u 
dt 

= ((7 + AY'Pu - (7 + A)' 13 (tu - V 2 (x, u)) , d t u) H0{NTd) 
= + Af-^n- ^ + A)-^ { 1 u-V 2 {x,u))\\ 2 H , [NTd) <0. 

Recalling that Ai < a(x) < A 2 , we have S N (u(t)) < S N (u ) < -^A 2 |w| 2 + 
A^II^Hl- for all t > 0. We can conclude that 

\Vu(x, t)\ 2 dx < [ -a(x)Vu(x, t) ■ Vu(x, t)dx 



[0,N] d J[0,N] d 2 



(33) 

< ^A^ 2 + 2N d \\V\\ L ^ N]d) . 



N d 



Thus, ||Vu(t)||^2 is bounded uniformly in t. 

Because Sn(u) is bounded below and ^<SW(u(i„)) < 0, there is a sequence t n — > 
00 such that ^SW(u(i n )) — > as n — > 00. The periodicity of V ensures that for any 
sequence of integers, {k n }, we have 5jv(u(f n ) +fe Ti ) = S]y(u(t n )). By selecting fc„ = 
— [N~ d Jj ^ u(x, t„)dx\ , we can replace the sequence u(t n ) by u(t n ) + k„ without 
affecting Sn- Therefore we assume, without loss of generality, that for each n, the 
average of u(x,t n ) lies in the range [0,1]. The gradients Vu(x,t n ) arc uniformly 
bounded in L 2 , so by Poincare's inequality we have that {u(t n )} is precompact, 
and there is a subsequence u(t nk ) that converges weakly in H 1 and strongly in L 2 . 
Denote the limit by u*. V 2 is Lipschitz, so ||V2(-,ti(t nfc )) — V 2 (-, u* )\\ L 2 — > as 
k — > 00. A is an I/ 2 -closed operator, hence Au* + V 2 (x, m* ) =0. □ 

Lemma 6.3. Ifu(x,t) solves wzi/i initial condition u(x, 0) = uo(x) G -B W) i/ien 
G S w /or a/Z f > 0. 



20 



T. BLASS, R. DE LA LLAVE, AND E. VALDINOCI 



Proof. As before, we denote the solution u with initial condition uq of ([3]) as 
u(x,t) = &tUo(x). We can conclude from Theorem 12.11 that if uq < vq then 
<&tUo(x) < ^tVo(x). For convenience we let and IZi denote the family of op- 
erators 

Cku(x) — u(x + k) and IZiu(x) = u(x) + I 

for each k £ 7L d and I £ Z. 

Let I £ Z and define Yi(x,t) = lZi&tUo{x) = u(x,t) + I. Recall the abbreviated 
notation 

Lu = -{ 1 + A) 1 - f3 u and X(u) = (7 + A)- p {^u - V 2 (x, u)), 

which allows us to write ((3]) as dtu = Lu + X(u), or dt&tUo = £$t«o + A(<£> t ito). 

Now, V 2 is periodic over the integers, so V 2 (x, u) = V 2 (x, TZiu). Also, (7 + A) Q c = 
"f a c for any constant c£l. Hence 

XiYi) = X(K&tua) =X(u(x,t)+l) = (7 + A)- (>yu + 7I - V 2 (x, u + I)) 
= 7 1_/3 « + (7 + A)-P{ 1U - V 2 (x, u)) = 7 1 - /3 ; + X($ t u ). 

Similarly, 

LYi = L(Ki$tU ) = L(u(x, t) + 1) = -(7 + A) l -^{u{x, t) + I) 
= - (7 + Af-^uix, t) = -j 1 -^ + L($ t u ). 

Therefore L$ t u + X($ t u ) = LY t + X(Yi) and we have d t Yi = d t § t uo = 
L$ t u + X($ t u ) = LYi + X{Y l ). So Y t solves © with initial condition Y t (x, 0) = 
Uq(x)+1. But d t ($t'R-i u o) = L(<f> t 1liUo)+X(<f> t 1ZiUo) with the same initial condition 
QqIZiUq = u + I. Thus, 

<5>tn lUo = Ui^tuo (34) 

by the uniqueness of solutions to ((3]) as shown in Proposition 13.61 

Now define Zk{x,t) — Ck&tuo = u(x + k,t). The periodicity of V 2 ensures that 
CkV 2 (x, u(x, t)) = V 2 (x + k, u{x + k, i)) = V 2 (x, u(x + k, t) so that X(Cku(x, t j) = 
CkX(u(x,t)). Clearly L(Cku(x,t)) = CkLu(x,t), so 

dtZk = dtu(x + k,t) = CkLu(x, t) + CkX(u(x, £)) 
= LC k u(x, t) + X(C k u(x, t)) = LZ k + X(Z k ). 

So Zk solves (|3]) with initial condition Zi(x, 0) = uq{x + k). But dt(&tCkUo) = 
L(&tCkUo) + X (<f>tCfeMo) with the same initial condition <f>oCfcito = uo(x + k). Thus, 

$ t C fc u = C fc $ t u (35) 

by the uniqueness of solutions to as shown in Proposition 13.61 

Now suppose uq £ B^ so that CuR-iUq < or > according to whether uj-k+l < 
or uj-k+l > O.Then the comparison principle from Thcorcm l2.1l vields $ t CkR-iUQ < 
or > according to whether uj-k + l<0 or uj-k + l>0. But equations (|3"4"|) and 
([35]) together give $ t Cfc7?./ = CkTZi&t, so CkR-i^tUo < or > according to whether 
uj ■ k + I < or lo ■ k + I > 0, hence $ t w £ B u . □ 

The set is closed under LMimits, so u* £ B^. This establishes Theorem 16. II 
in the case u> £ -hZ/ 1 . 
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6.2. Irrational Frequencies. For the case of irrational frequency, let ui £ M. d \Q d 
and let (uj n ) be a sequence such that w„ £ \l* d for each n and u> n — > uj. Let 
div(a(a;)Vu„) = V2{x,u n ), where u n £ uj n ■ x + if 1 (7VT d ), A < a(x) < A, and let 
Br denote a ball of radius R centered at some point in the domain Q such that the 
concentric ball B±r is also in the domain f2. Then Theorem 8.22 of |GT01j gives 
for each n that 

osc Bh (m„) < (1 - C~ 1 )osc BlR {u n ) + k(R) 

where k = X^ 1 R 2 ^ 1 ~ d / q ^\\V2{x, u n )\\Li/ 2 (n) f° r 1 > d, which is bounded by the 
volume of ft times ||F 2 |U~, and C = C(d,A/X,q) > 1 (pag es 200-2 01 of |GT01| l 
Following Moser's methods from the proof of Theorem 2.2 of Mos86 , we can obtain 
a local result on the cube Q = {x £ R d : \x 3 ■ | < 1 /2}, and then consider translations. 
Using a variant of the above inequality we can say 

oscq(u„) < 6>osc 4 q(u„) + ci 

where 9 £ (0, 1) depends on A, A, d, since we can fix q = d + 1, and ci depends 
on A, d, and || V~2{x, u n )\\ Lq /2^ n y Here we will take fl = 5Q = {\xj\ < 5/2} so that 
\\V 2 {x,u n )\\ Lq/ 2 in) < 5 d ||(y 2 )( d+1 )/ 2 || L oo, independent of n. 

Since u„(x + k) — u n {x) — / has the same sign as uj n ■ k — I, for any x £ M. d and 
any k £ Z d , I £ Z, we can conclude \u n {x + k) — u n (x) — uj„ ■ k\ < 1 for any such x 
and k. For an arbitrary y £ M. d select k £ Z d such that y — k £ Q, then we have 

\u n (x + y) - u n (x) - y ■ u n \ < \u n (x + y) - u n (x + k)\ + 

+ \u n (x + k) - u„(x) - uj n ■ k\ + \uj„ ■ k - u n -y\ 
< osc x+k +Q(u n ) + 1 + \u) n \ 

Then an argument parallel to the one on page 240 of |Mos86j shows that osc4q(w„) < 
oscq(u„) + 22 J= i(l + 2|(w„)j|) so that with the result from |GT01| we have 
oscq(u„) < C2-JI + \oJ n \ 2 - This estimate holds in translated cubes x + Q, so if 
we assume that |w n | < c\ui\, then the \u n {x + y) — u n (x) — uj n ■ y\ are bounded in 
R d by 031/1 + \oj\ 2 , independently of n. 

This bound is crucial, because it allows the application of Theorem 5.2, page 277 
of [LU68I. which gives u n £ u) n ■ x + C l e (Assuming V £ C 2,e ), and bounds 

|Vu„| < 71 

with 7! depending on the ellipticity constants and cu n (in fact, monotone in w„), so 
bounded for our sequence of (w„). So, convergence of tu n — > uj and the Arzela-Ascoli 
Theorem imply the existence of a subsequence u nk — * u* in C;° oc . Then, div(a(a;)V) 
is closed under C ; ° oc limits, and we have 

Au* u + V 2 (x, u*J = -div(a(x)V<) + V 2 (x, u*J = 0. 

This completes the proof of Theorem 16.11 □ 

Remark. If wc replace the operator A from Theorem 16.11 with A a for a > 0, 
then the same method of proof for Theorem 16.11 will allow us to find solutions of 
A a u + V2{x,u) = 0, with u £ provided that uj is rational. However, to find 
solutions for irrational w, we require a oscillation lemma for the critical points with 
rational frequencies. With this, one may be able to find solutions with irrational 
frequencies using methods from [CSM05| . The work of |CS09| establishes a Harnack 
inequality in the autonomous case, and can perhaps be extended to our case as well. 
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For this method to work, one would need both the oscillation lemma of De 
Giorgi-Moser type as well as C e estimates for the fractional power of an elliptic 
operator. 
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